Abrikosov Lattice Solutions of the Ginzburg- Landau Equations 

T. Tzaneteas*t and I. M. Sigal « 

Dept. of Mathematics, Univ. of Toronto, Toronto, Canada, M5S 2E4 

October 31, 2009 

Abstract 

Building on the earlier work of Odeh, Barany, Golubitsky, Turski and Lasher we give a proof of existence of 
Abrikosov vortex lattices in the Ginzburg-Landau model of superconductivity. 
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1 Introduction 

1.1 The Ginzburg-Landau Model. The Ginzburg-Landau model of superconductivity describes a superconductor 
contained in f2 C R", n = 2 or 3, in terms of a complex order parameter ip : f2 — * C, and a magnetic potential 
A : — * Kxl The key physical quantities for the model are 

• the density of superconducting pairs of electrons, n s := \ip\ 2 ', 

• the magnetic field, B := curl A; 

• and the current density, J :— Im^Y^). 

In the case n = 2, curl^4 := ^A 2 - — ^A 1 - is a scalar- valued function. The covariant derivative Va is defined to be 

' OXl 0X2 

V — iA. The Ginzburg-Landau theory specifies that a state (ip, A), in the absence of an external magnetic field, has 
energy 

£ n (lM):= / |V^| 2 + |curlA| 2 + ^(l-M 2 ) 2 , (1) 
Jn 1 

where k is a positive constant that depends on the material properties of the superconductor. 

It follows from the Sobolev inequalities that for bounded open sets ft, £q is well-defined and C°° as a functional on 
the Sobolev space H . The critical points of this functional must satisfy the well-known Ginzburg-Landau equations 
inside ft: 

A A V = ^(iVf (2a) 



curl*curlA = Im('0V4'0)- (2b) 

Here — —V£^a, V| and curl* are the adjoints of ^4 and curl. Explicitly, V^F = — div F + iA ■ F, and 
curl* F = curlF for n = 3 and curl* f = f-J^-, — -J^-) for n = 2. 

J v OX2 ' OXl J 

There are two immediate solutions to the Ginzburg-Landau equations that are homogeneous in tp. These are the 
perfect superconductor solution where "0s = 1 and As = 0, and the normal (or non-superconducting) solution where 
tpN = and An is such that curl A at =: b is constant. 
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1 The Ginzburg-Landau theory is reviewed in every book on superconductivity. For reviews of rigorous results see the papers 1101 1111 
and the books [231 fl5l [Tfl 122] 
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It is well-known that there exists a critical value k c (in the units used here, k c = 1/V2), that separates supercon- 
ductors into two classes with different properties: Type I superconductors, which have n < k c and exhibit first-order 
phase transitions from the non-superconducting state to the superconducting state, and Type II superconductors, 
which have K > k c and exhibit second-order phase transitions and the formation of vortex lattices. Existence of the 
vortex lattice solutions is the subject of the present paper. 

1.2 Results. In 1957, Abrikosov pQ discovered solutions of whose physical characteristics n s , B, and J 
are periodic with respect to a two-dimensional lattice, while independent of the third dimension, and which have a 
single flux per lattice celH. (In what follows we call such solutions, with n s and B non-constant, lattice solutions, 
or, if a lattice C is fixed, C-lattice solutions. In physics literature they are called variously mixed states, Abrikosov 
mixed states, Abrikosov vortex states.) Due to an error of calculation he concluded that the lattice which gives the 
minimum average energy per lattice cell is the square lattice. Abrikosov's error was corrected by Kleiner, Roth, and 
Autler [18] , who showed that it is in fact the triangular lattice which minimizes the energy. 

Since then these Abrikosov lattice solutions have been studied in numerous experimental and theoretical works. 
Of more mathematical studies, we mention the articles of Eilenberger [14] and Lasher [19] . 

The rigorous investigation of Abrikosov solutions began soon after their discovery. Odeh [21] proved the existence 
of non-trivial minimizers and obtained a result concerning the bifurcation of solutions at the critical field strength. 
Barany, Golubitsky, and Tursky [5] investigated this bifurcation for certain lattices using equivariant bifurcation 
theory, and Takac [24] has adapted these results to study the zeros of the bifurcating solutions. 

Except for a variational result of [ST] (see also [13]), work done by both physicists and mathematicians has 
followed the general strategy of [JJ . 

In this paper we combine and extend the previous technique to give a self-contained proof of the existence of 
Abrikosov lattice solutions. To formulate our results we mention that lattices £cR 2 are characterized by the area 
\Qc\ of the lattice cell ftc and the shape r, given by the ratio of basis vectors identified as complex numbers (for 
details see Section [3]). We will prove the following results, whose precise formulation will be given below (Theorem 

El). 

Theorem 1. Let C be a lattice with \ \flc \ — lf| *C 1. 

(I) If |fic| > =r, there exists an C-lattice solution. If \0lc\ < =5, then there is no C-lattice solution in a neigh- 
bourhood of the branch of normal solutions. 

(II) The above solution is close to the branch of normal solutions and is unique, up to symmetry, in a neighbourhood 
of this branch. 

(III) The solutions above are real analytic in \£lc\ in a neighbourhood o/=$. 

(IV) The lattice shape for which the average energy per lattice cell is minimized approaches the triangular lattice as 
1^1- 

Remark 2. 

(a) [3TJ[TJ] showed that for all |Q,c| > ^ there exists a global minimizer of £q c . 

(b) [HJ[S] proved results related to our solutions in (I). 

(c) [H] proved partial results on (IV). 

Among related results, a relation of the Ginzburg-Landau minimization problem, for a fixed, finite domain and for 
increasing Ginzburg-Landau parameter k and external magnetic field, to the Abrikosov lattice variational problem 
was obtained in [3[[5]. [12] (see also [13] ) have found boundaries between superconducting, normal and mixed phases. 

All the rigorous results above deal with Abrikosov lattices with one quantum of magnetic flux per lattice cell, 
partial results for higher magnetic fluxes were proven in [5] 3] . This problem will be addressed in our subsequent 
paper. 

Acknowledgements 
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I.M.S.'s stay at the IAS, Princeton. 



2 Such solutions correspond cylindrical samples. In 2003, Abrikosov received the Nobel Prize for this discovery 



AbrikosovLattices, October 31, 2009 



3 



2 Properties of the Ginzburg-Landau Equations 

2.1 Symmetries. The Ginzburg-Landau equations exhibit a number of symmetries, that is, transformations which 
map solutions to solutions. The most important of these symmetries is the gauge symmetry, defined for any sufficiently 
regular function r\ : — > R, which maps (tp, A) i— > (T^ip, T^A), where 

T„V = e"V, T n A = A + Wi]. (3) 

There are also the translation symmetry, defined for each t 6 R 2 , which maps (ip, A) i— > (T t ip,T t A), where 

T t ip{x) := i){x + 1), T t A(x) := A(x + t), (4) 

and rotation and reflection symmetry, defined for each R £ 0(2) (the set of orthogonal 2x2 matrices), which maps 
0>,A) ^ (Tr^TrA), where 

T R ip(x) :=i()(Rx), TrA{x):=R- 1 A{Rx). (5) 

2.2 Flux Quantization. One can show that under certain boundary conditions (e.g., 'gauge-periodic', see 
below, or if Q, = R 2 and £q < oo) the magnetic flux through il is quantized. 

3 Lattice States 

Our focus in this paper is on states (ip,A) defined on all of R 2 , but whose physical properties, the density of 
superconducting pairs of electrons, n s := \ip\ 2 , the magnetic field, B := curl A, and the current density, J := 
Im(^V4^), are doubly-periodic with respect to some lattice C. We call such states C— lattice states. 

One can show that a state (ip, A) £ Hl oc (M. 2 ; C) x Hl oc (R 2 ; R 2 ) is a £-lattice state if and only if translation by an 
element of the lattice results in a gauge transformation of the state, that is, for each t G C, there exists a function 

9t e ffL( R2 ; R ) such that 

f(x + t) = e l9t{x) iP(x) and A(x + t)= A(x) + Vg t (x) 

almost everywhere. 

It is clear that the gauge, translation, and rotation symmetries of the Ginzburg-Landau equations map lattice 
states to lattice states. In the case of the gauge and translation symmetries, the lattice with respect to which the 
solution is periodic does not change, whereas with the rotation symmetry, the lattice is rotated as well. It is a 
simple calculation to verify that the magnetic flux per cell of solutions is also preserved under the action of these 
symmetries. 

Note that (ip, A) is defined by its restriction to a single cell and can be reconstructed from this restriction by 
lattice translations. 

3.1 Flux quantization. The important property of lattice states is that the magnetic flux through a lattice cell 
is quantized: 

$(A) := / curl A = 2vm (6) 
Jn 

for some integer n. Here is any fundamental cell of the lattice. Indeed, if > on the boundary of the cell, we 
can write i/j = \ip\e and < 9 < 2tt. The periodicity of n s and J ensure the periodicity of V(9 — A and therefore by 
Green's theorem, J n curl A = <f dn A = <f gn W8 and this function is equal to 27m since ip is single- valued. 

We let b be the average magnetic flux per lattice cell, b := j^$(A). Equation |6]) then imposes a condition on 
the area of a cell, namely, 

Due to the physical interpretation of b as being related to the applied magnetic field, from now on we use b = as 
a parameter of our problem. We note that due to the reflection symmetry of the problem we can assume that b > 0. 

3.2 Lattice Shape. In order to define the shape of a lattice, we identify x € R 2 with z — x\ + 1x2 <E C, and 
view £ as a subset of C. It is a well-known fact (see [6]) any lattice £CC can be given a basis r,r' such that the 
ratio t — ^- satisfies the inequalities: 

(i) M > 1. 
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(ii) Imr > 0. 

(iii) -\ < Rer < ±, and Rcr > if |t| = 1. 

Although the basis is not unique, the value of r is, and we will use that as a measure of the shape of the lattice. 

Using the rotation symmetry we can assume that if C has has as a basis {rei,rr}, where r is a positive real 
number and e\ = (1,0). 

3.3 Fixing the Gauge. The gauge symmetry allows one to fix solutions to be of a desired form. We will use 
the following preposition, first used by [H] and proved in [53] . We provide an alternate proof in Appendix [Cl 

Proposition 3. Let (tp,A) be an C-lattice state, and let b be the average magnetic flux per cell. Then there is a C- 
lattice state (4>, ^4o+ a ) that is gauge-equivalent to a translation of (ip, A), where Aq(x) = \x^- (where x 1 - = (— x%, x\) ), 
and <j) and a satisfy the following conditions. 

(i) a is doubly periodic with respect to C: a(x + 1) = a(x) for all t € C. 

(ii) a has mean zero: J fl a = 0. 

(iii) a is divergence-free: diva = 0. 

(iv) (j>(x + t) = e^ tAx (f)(x), where t A x — tix% — tix\, for t = re\,rr. 

3.3 Lattice Energy. Lattice states clearly have infinite total energy, so we will instead consider the average 
energy per cell, defined by 



Here, fi is a primitive cell of the lattice with respect to which (ip, A) is a lattice state and \ is its Lebesgue measure. 
We seek minimizers of this functional under the condition that the average magnetic flux per lattice cell is fixed: 

In terms of the minimization problem, we see that the perfect superconductor is a solution only when <i>(A) is 
fixed to be = 0. On the other hand, there is a normal solution, (ipx = 0, An, curl An = constant), for any 

condition on Q(A). 

We define the energy of the lattice with the flux n per cell as 



where the infimum is taken over all smooth /^lattice states satisfying (i) through (iv) of Proposition [3] 

3.4 Result. Precise Formulation. The following theorem gives the precise formulation of Theorem [T] from 
the introduction. 

Theorem 4. Let n = 1. 

(I) For every b sufficiently close to but less than the critical value b c — k 2 , there exists an C— lattice solution of 
the Ginzburg- Landau equations with one quantum of flux per cell and with average magnetic flux per cell equal 
to b. 

(II) This solution is unique, up to the symmetries, in a neighbourhood of the normal solution. 

(III) The family of these solutions is real analytic in b in a neighbourhood of b c . 

(IV) If k 2 > 1/2, then the global minimizer Lb of the average energy per cell, £i(C), approaches the Ctrianguiar as 
b — > b c in the sense that the shape Tf, approaches T tr i angu i ar — e lTr / 3 in C. 

The rest of this papers is devoted to the proof of this theorem. 



Ety,A) :=—£n(i>,A). 



(7) 



£ n {C) :=iafE(il>,A), 



(8) 
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4 Rescaling 

In this section we rescale (ip,A) to eliminate the dependence of the size of the lattice on b. Our space will then 
depend only on the number of quanta of flux and the shape of the lattice. 

Suppose, then, that we have a £-lattice state (ip,A), where C has shape r. Now let b be the average magnetic 
flux per cell of the state and n the quanta of flux per cell. From the quantization of the flux, we know that 

b 



|f2| r 2 Imr' 

We set a := (f) 2 ■ The last two relations give a — (^r) 2 r. We now define the rescaling (ip, A) to be 

(tp(x),A(x)) := (cnp(ax),aA(ax)). 
Let £ T be the lattice spanned by r T and r r r, with fl T being a primitive cell of that lattice. Here 

_ T _ ( 27T 



We note that \il T \ = 2nn. We summarize the effects of the rescaling above: 

(i) (ip,A) is a £ T -lattice state. 

(ii) E(ijj,A) = £x(i>,A), where A = ^ and 

ix{ ^ A) = ^>?L (l v ^l 2 + l curM l 2 + y(^l 2 -^) 2 ) **■ ( 10 ) 

(iii) ip and A solve the Ginzburg-Landau equations if and only if and A solve 

(-A A - \)t/> = -« 2 |V|V, ( lla ) 

curl* curl A = lm{{/;\7jiip} (lib) 

2 

for A = The latter equations are valid on fi T with the boundary conditions given in the next statement. 

(iv) If (ip, A) is of the form described in Proposition then 

A = Aq + a, where A%(x) := ^ x± : 

where x 1 - = (— X2, Xi), and ip and a satisfy 

(a) a is double periodic with respect to C T , 

(b) /„ T a = 0, 

(c) div a = 0, 

(d) ^(a; + 1) = e^ tAx ip(x) for t = r r , r T r. 

In w/iat follows we drop the hat from ip, A, and £\. 

We are now state our problem in terms of the fields ip and a. We define the Hilbert space % n (r) to be the 
closure under the L 2 -norm of the space of all smooth ip on fl T satisfying the quasiperiodic boundary condition (d) 
in part ([Iv]) above. M t n {r) is then the space of all ip S S£ n (r) whose (weak) partial derivatives up to order 2 are 
square-integrable. 

Similarly, we define the Hilbert space S£(r) to be the closure of the space of all smooth a on fl T that satisfy 
periodic boundary conditions, have mean zero, and are divergence free, and J4?{t) is then the subspace of _S? (r) 
consisting of those elements whose partial derivatives up to order 2 are square-integrable. 

Our problem then is, for each n — 1,2,..., find (ip, a) £ <^ n ( T ) x ^( T ) so that (ip, Aq + a) solves the rescaled 
Ginzburg-Landau equations (fTTj) . and among these find the one that minimizes the average energy £\. 
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5 Reduction to Finite-dimensional Problem 

In this section we reduce the problem of solving Eqns to a finite dimensional problem. We address the latter in 
the next section. Substituting A = Aq + a, we rewrite (|IIj) as 

(L n - X)ip + 2ia • Va^ + \a\ 2 ip + k 2 \^\ 2 i}j = 0, (12a) 

(M+|^| 2 )a-Im{^V4^}-0, (12b) 

where 

L n := -A A? and M := curl* curl . (13) 

The operators L n and M are elementary and well studied. Their properties that will be used below are summarized 
in the following theorems, whose proofs may be found in Appendix IB] 

Theorem 5. L n is a self-adjoint operator on J^ n (r) with spectrum u{L n ) — {{2k + l)n : k — 0,1,2,...} and 
dime null(i™ — n) = n. 

Theorem 6. M is a strictly positive operator on M 7 (t) with discrete spectrum. 

We first solve the second equation (|12b[) for a in terms of ip, using the fact that M is a strictly positive operator, 
and therefore M + \ip\ 2 is invertible. We have a = a(ip), where 

a(V0 = (M + IVI 2 )- 1 Im(^VA^). (14) 

We collect the elementary properties of the map a in the following preposition, where we identify with a real Banach 
space using ip «-» ip := (Re ?/>,Im tp). 

Proposition 7. The unique solution, a(ijj), of (I12b|) maps J^ n (r) to Jf?(r) and has the following properties: 

(a) a(-) is analytic as a map between real Banach spaces. 

(b) o(0) = 0. 

(c) For any a£l, a(e' iQ V) = ^(ip)- 

Proof. The only statement that does not follow immediately from the definition of a is (a). It is clear that Im^V*™ VO 
is real-analytic as it is a polynomial in ip and Vi/ 1 , and their complex conjugates. We also note that {M — z) _1 is 
complex-analytic in z on the resolvent set of M, and therefore, (M + IV'I 2 ) -1 is analytic, (a) now follows. □ 

Now we substitute the expression (fl4|) for a into (|12a|) to get a single equation F(\,ip) — 0, where the map 
F:lxJf„(T)^ Jz?„(t) is defined as 

F(A, V) = (L n - A)V + 2ia{i>) • Va»V + + « 2 |^|V (15) 

The following proposition lists some properties of F. 
Proposition 8. 

(a) F is analytic as a map between real Banach spaces, 

(b) for all X, F(A,0) = 0, 

(c) for all X, L> v ,F(A,0) = L n - A, 

fdj for all «el, F(A, e M V) = e ia F(X, ip). 
(e) for all ip, (ip,F(X,t/})) G M. 
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Proof. The first property follows from the definition of F and the corresponding analyticity of a(ip). (b) through (d) 
are straightforward calculations. For (e), we calculate that 

(iP,F(X,iP)) = (iP,(L n -X)iP)+2i f ijja(iP)-ViP + 2 [ (c#) • A )\iP\ 2 + f |aW>)| 2 |Vf + k 2 / IV>| 4 - 

JSJT 7[Jt Jq t Jq t 

The final three terms are clearly real and so is the first because L n — A is sclf-adjoint. For the second term we 
calculate the complex conjugate and see that 



2i / tpa(ip) ■ Vip = -2i / ipa(ip) ■ Vtp = 2i / (V^ • ot(ip))ip, 
Jsi T Jn-r Jq. t 

where we have integrated by parts and used the fact that the boundary terms vanish due to the periodicity of the 
integrand and that div a(ip) = 0. Thus this term is also real and (e) is established. □ 

Now we reduce the equation F(X,tp) = to an equation on the finite-dimensional subspace null(L" — n). To 
this end we use the standard method of Lyapunov-Schmidt reduction. Let X := J^f n {r) and Y :— J£ N (r) and let 
K = null(L" — n). We let P be the Riesz projection onto K, that is, 

P~--LL L n -z)- 1 dz, (16) 

27TI J 7 

where 7 C C is a contour around that contains no other points of the spectrum of L. This is possible since is an 
isolated eigenvalue of L. P is a bounded, orthogonal projection, and if we let Z := nullP, then Y = K © Z. We also 
let Q :— I — P, and so Q is a projection onto Z. 

The equation -F(A, if>) = is therefore equivalent to the pair of equations 

PF(X, Pi/j + Qip) = 0, (17) 
QF{\,Pxj; + Q^) = 0. (18) 



We will now solve (|18|) for w = Qip in terms of A and v = Pip. To do this, we introduce the map G : M.X K x Z — > Z 
to be G(A, v, w) := QF(X, v + w). Applying the Implicit Function Theorem to G, we obtain a real-analytic function 
w : R x K — > Z , defined on a neighbourhood of (n, 0), such that w = w(X, v) is a unique solution to G(A, v, w) = 0, for 
(A, v) in that neighbourhood. We substitute this function into (|17[) and see that the latter equation in a neighbourhood 
of (n, 0) is equivalent to the equations 

tp = v + W (X, v) (19) 

and 



7 (A,w) := PF(X,v + w{X,vj) = (20) 

(the bifurcation equation). Note that 7 : R x X — > C. We have shown that in a neighbourhood of (n, 0) in R x X, 
(A, ip) solves F(X, tp) = if and only if (A, v), with v — Pip, solves (I2TJ1) . 
Finally we note that 7 inherits the symmetry of the original equation: 

Lemma 9. For every a£l, 7(A, e m t;) = e za j(X,v). 

Proof. We first check that w(X, e la v) — e za w(X,v). We note that by definition of w, G(X,e la v,w(X,e la v)) — 0, but 
by the symmetry of F, we also have G(A, e la v, e la w(X, v)) = e la G(X,v,w(X,v)) — 0. The uniqueness of w then 
implies that w(X, e la v) = e la w(X, v). We can now verify that 

7 (A, e la v) = PF(X, e la v + w(X, e ia v)) = e la PF(X, v + w(X, v))) = e M 7 (A, v). □ 

Solving the bifurcation equation P0|) is a subtle problem unless n = 1. The latter case is tackled in the next 
section. 
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6 Bifurcation Theorem, n = 1 

In this section we look at the case n = 1, and look at solutions near the trivial solution. For convenience we drop 
the index n = 1 from the notation. We will see that as b = ^- decreases past the critical value b = k , a branch of 
non-trivial solutions bifurcates from the trivial solution. More precisely, we have the following result. 

Theorem 10. For every t there exists a branch, (X s ,ip s , A s ), s G C with \s\ 2 < e for some e > 0, of nontrivial 
solutions of the rescaled Ginzburg- Landau equations unique (apart from the trivial solution (\,0,A<f)) in a 

sufficiently small neighbourhood of (1,0, Aq) in R x J4?(t) x ,^(t), and s.t. 

'\ s = l+g x (\s\ 2 ), 
ip s = sip + sg^ds] 2 ), 
^A s = A Q +g A (\s\ 2 ), 

where (L — l)^o = 0, g^ is orthogonal to null(L — 1), g\ : [0,e) — > R, g^ : [0, e) — > ,ffl(r), and gA ■ [0,e) — > ^(r) 
are real-analytic functions such that g\(0) = 0, g,/,(0) = 0, 5a (0) = and 3^(0) > 0. Moreover, 



Proof. The proof of this theorem is a slight modification of a standard result from the bifurcation theory. It can be 
found in Appendix VK[ Theorem [15l whose hypotheses are satisfied by F as shown above (see also [21] [8]). The latter 
theorem gives us a neighbourhood of (1,0) in R x ^(r) such that the only non-trivial solutions are given by 

fA s = l + . 9A (| s | 2 ), 
|Vs = sipo + sg^(\s\ 2 ). 

Recall a(ip) defined in (| 14|) . We now define c/a(s) — a(ip s ), which is real-analytic and satisfies §a(— t) = a (~V-'t) — 
gA(t), and therefore is really a function of t 2 , gA(t 2 )- Hence A s = Aq + gA(\s\ 2 ). 

Finally, to prove (|2ip we multiply the equation ^(A, if) = scalarly by ipo and use that L is self-adjoint and 
(L — 1)i/jq = to obtain 

(V>o, (A - = 2i(^ , a(V) ' Va V) + <^o, |aty>)|V) + « 2 <V>o, MV>- 
Let a\ : = g' A (0). Substituting here the expansions obtained in the first part of the theorem, we find 

g' x (0)\\M 2 = ^ f $oai • V Aa ^ + K 2 ( |^o| 4 . (22) 

In order to simplify this expression we first note that by differentiating (| 12b|) w.r. to \s\ 2 at s = 0, we obtain 
curl* curlai = Im(?/'oV4oV'o) ■ 

Now for the first term on the r.h.s. of ([22]) , taking the imaginary part of (|22]) we see that Re(/ nT t/>o«i " ^loV'o) = 
and therefore 



2i I V ai • VaqV'o = -2 / a x • Im(V>oY4 o '0o) ^ 
= —2 / ai • curl* curlai dx 



-2 / (curlai) 2 dx. 



Here in the second step we integrated by parts. Next we show that 

-icurl* |^o| 2 = M^oVaoVo)- ( 23 ) 
Using the notations of Appendix 151 (with n = 1), we have that L-ipo = 0: 

1 i 
dx^o + id X2 if> + -xi-0o + ^x 2 ipo = 
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Multiplying this relation by ipo and subtracting and adding the complex conjugate of the result, we obtain the two 
relations 



This means that 
Im(^oV4 ^o) 



('4 , odx 2 ipo ~ i>od X2 ipo) - § ari |Vo| 2 



%_ 

2 

-5(^0^x2^0 + ipod X2 ^o) 



i9 X2 \M 2 







which gives (|23p. 

(j2"5)) implies that curlai = — llV'ol 2 + C for some constant C. This C can be determined using the fact that, since 
A has mean zero, cti does as well, and this gives C = ^ J QT IV'ol 2 , which establishes 

curla 1 = -i|^ | 2 + ^y'jV'o| 2 - (24) 
Using this equation we finish the calculation above: 

2iJ^^oai-S/ Ao ipodx = ~J^ |^ | 4 Gb+^(y |^o| 2 ^ . (25) 

Substituting this expression into (|22"|) and rearranging terms we arrive at (f2"l"|) . And that completes the proof of 
Theorem OH □ 

Theorem rrU] implies (I) - (III) of Theorem S □ 
Finally, we mention 

Lemma 11. Recall that Imr > 0. Let (A s , ip s , A 3 ) be the solution branch constructed above and let m T — 

/ 1 T \ _ _ 

(Vlmr) -1 (oi )' ^ erl (^s,ifis, A s ), where the functions (ip s ,A s ) are defined on a T-independent square 

lattice and are given by 

{$ s (x) = ip s (m T x), 
\A s (x) = M$A a (m T x), 

depend R- analytically on t . 

We sketch the proof of this lemma. The transformation above maps functions on a lattice of the shape r into 
functions on a T-independent square lattice, but leads to a slightly more complicated expression for the Ginzburg- 
Landau equations. Namely, let U T ip(x) :— tp(m T x) and V T a(x) := m T a(m T x). Applying U T and V T to the equations 
(j!2a[) and (|12b|) , we conclude that (ip s , A s ) satisfy the equations 

(£? - X)tp + 2i(m t T )- 1 a • (mt)- 1 V A »V + + k 2 |V>|V = 0, (27a) 



(M r + |V| 2 )a + ^ T Q (V) =0, (27b) 

where 

L" := -UrAA^U' 1 and M T := U r curl* curlF,. -1 . (28) 

Here we used that V t Aq = A£ and U T Vij) = (m^) _1 r/ r '0- (The latter relation is a straightforward computation 
and the former one follows from the facts that for any matrix m, (mi) 1 = (det m)(m t )~ 1 x ± , and that in our 
case, detm T = 1.) Note that the gauge in the periodicity condition will still depend on Imr. These complications, 
however, are inessential and the same techniques as above can be applied in this case. The important point here is 
to observe that the function ipg, constructed in Appendix B, the function w(X, sipo), where w(\, v) is the solution 
of (fT8|) , and the bifurcation equation (f20|) depend on r real-analytically. We leave the details of the proof to the 
interested reader. 
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7 Abrikosov Function 

In this section, we continue with the case n = 1. We prove the Abrikosov relation between the energy per cell and 
the Abrikosov function, (3{t), which is defined by 

m := TT^W' (29) 

where ipQ is a non-zero element in the nullspace of L n — 1 acting on J^(r). Since the nullspace is a one-dimensional 
complex subspace, (3 is well-defined. 

Recall that b = Since the function <?a(|s| 2 ) given in Theorem ITD1 obeys g\(0) = and <? A (0) ^ 0, the function 
b s = k 2 (1 + ,<7a (I s ! 2 )) 1 =: ^ 2 + .°f>(l s | 2 ) can be inverted to obtain \s\ = s(b). Absorbing s = ]fy into we can define 
the family (^J^s , -^g(&) j ) °^ £ T_ periodic solutions of the Ginzburg-Landau equations parameterized by average 
magnetic flux 6 and their energy 

E b (r) := £ K 2 /b (i}j T s{b)l A T s(b) ). 

Clearly, , ^s(f>) > ^s(b) are analytic in 6. We note the relation between the new perturbation parameter := n 2 — b 
and the bifurcation parameter |s| 2 : 

M=^ 2 =» 2 H 2 + 0(| S r). (30) 

The relation between the Abrikosov function and the energy of the Abrikosov lattice solutions is as follows. 
Theorem 12. In the case n > the minimizers, t\>, of t h- > 2?(,(t) are related to the minimizer, t*, of /3(t), as 



Tb — r* = 0(/W 2 ), In particular, Tb — > r* as 6 — > k 2 . 
Proof. We first show that the theorem is a consequence of the following proposition, which is proved below. 
Proposition 13. We /iave 

*M = ^ + - W", + (£ - 1 + m ;_, Wt) ) + <V). (3D 

To prove the theorem we note that Eb(r) is of the form Eb(r) = eo + ei/i + e2(r)/i 2 + (3(/i 3 ). The first two 
terms are constant in r, so we consider Eb(r) — ea(r) + 0(/i). tj, is also the minimizer of t h £*(t) and t*, 
of ea(r). We have the expansions £&(r*) - Eb(n) = |-E b (t;,)(t* - r b ) 2 + 0((r* - n) 3 ) and E b (T*) - E b {n) = 
— \e 2 (Tb)(T*~Tb) 2 + 0((T* —Tb) 3 )+0([i), which imply the desired result. That concludes the proof of the theorem. □ 

Prof of Proposition \7M Recall that [i := k 2 — b. Using the real-analyticity of the function gb, gtfj, and g T A , we can 
express := n 2 /b, ^ T {y) := ipL^ and A T (p) := A T s , h s as 

A(/i) = l + iM + 0(M 2 ) (32) 
fa 

^(P) = M V Vo T + |U 3/ Vi + 0(M 5/2 ) (33) 

A T ( f x)=A Q +^a T 1 +0(ii 2 ). (34) 

We will first show that 

E b (r) = y + k '- + ^ (l - ^ j nT l^| 2 ) M 2 + 0{^). (35) 
Multiplying Qllap scalarly by ip and integrating by parts gives 

\V A i> T \ 2 =K 2 [ (A|r| 2 -K 2 ^ T | 4 ). 
Substituting this into the expression for the energy we find 
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Using the expansion above gives 

E b (r) = y + k4 - 2k ^ + £ ( X - Y J l^o I 4 + J I curlal | 2 ) ,, 2 + 0(/* 3 ), (37) 
where we have used the fact that curlAo = 1. As we proved above in (|24|) . we can show that 

curla[ = ~\\% I 2 + ^j n y \ 2 - 
Substituting this expression into (f37|) . we obtain 

E b (r) = y + - 4 - 2*1* + ^ (l - (k 2 - ±) £ |^| 4 - ± (J^ |VoT) ^ ^ + 0(/i 3 ). (38) 

Now, if we differentiate (|12a,|) twice w.r. to y}l 2 at /i = to obtain 

(L - 1}4>1 = \% ~ « 2 |#| a # - 2za[ • V Ao ^. 
Now using the fact that L — 1 is self-adjoint and that (L — l)"0o = 0, we find 
= [ V>q(£ — l)ip{ dx 
= / ^{^l-^W^^-^l-VAj^dx 
= \ [ | I 2 dx - n 2 f |^| 4 dx - [ i>l(2ia\ ■ V Ao V>J) dx. 
An analogous calculation to the one in the proof of Theorem [10] then gives 

0=4/ \Vo\ 2 dx-K 2 [ \r \ A dx+\[ \^dx--}-([ \Vo\ 2 dx 

This relation gives ([33)) from Q38p. but also by dividing by (/ |"0ol 2 ) 2 an d rearranging we then obtain ([3"Tj) . □ 

The following result was discovered numerically in the physics literature and proven in [5] using earlier result of 
[20J: 

Theorem 14. The function /3(r) has exactly two critical points, r = e 1 *^ 3 and r = e m l 2 . The first is minimum, 
whereas the second is a maximum. 

Theorems [TOl [T2l [T4l after rescaling to the original variables, imply Theorem [4] which, as was mentioned above, 
a precise restatement of Theorem [T] of Introduction. 

A Bifurcation with Symmetry 

In this appendix we present a variant of a standard result in Bifurcation Theory. 

Theorem 15. Let X and Y be complex Hilbert spaces, with X a dense subset ofY, and consider a map F : Rx X ^ Y 
that is analytic as a map between real Banach spaces. Suppose that for some Xq 6 R, the following conditions are 
satisfied: 

1. F(X, 0) = for all X e K, 

2. D^,F(Xq,0) is self-adjoint and has an isolated eigenvalue at of (geometric) multiplicity 1, 

3. For non-zero v G null D^F(X , 0), (v,D\^F(X ,0)v) ^ 0, 

4. For all aeR, F(X,e ia tp) = e ia F(X, ip) . 

5. For all ipeX, (ip, F(X, ip)) € R. 
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Then (Ao,0) is a bifurcation point of the equation F(X,ip) — 0. In fact, there is a family of non-trivial solutions, 
(A, ip), unique in a neighbourhood of (Ao, 0) in R x X , and this family has the form 

Ja = <MN 2 ), 

\V> = sv + s(j)^{\s\ 2 ), 

for s G C with \s\ < e, for some e > 0. Here v G nullD^i^Ao, 0), and 4>\ : [0, e) — > R and 4>iji '■ [0, e) — > X are unique 
real-analytic functions, such that (f>\(0) — Xq, 4>^(0) = 0. 

Proof. The analysis of Section [S] reduces the problem to the one of solving the bifurcation equation PO)) . Since the 
projection P, defined there, is rank one and self-adjoint, we have 

p i' = Tr-^(v,ip)v, with v G nullL>^F(A ,0). (39) 

IMr 

We can therefore view the function 7 in the bifurcation equation (|20p as a map 7 : R x C — > C, where 

7(A, s) = (v, F(X, svq + w(X, sv)). 

We now look for non-trivial solutions of this equation, by using the Implicit Function Theorem to solve for A in 
terms of s. Note that if 7(A,i) = 0, then j(X,e la t) — for all a, and conversely, if 7(A, s) = 0, then 7(A, \s\) = 0. 
So we need only to find solutions of 7(A, t) — for (el. We now show that 7(A, t) £ R. Since the projection Q is 
self-adjoint, and since Qw(X,v) — w(X,v) we have 

(w(X, tv),F(X, tv + w(X, tv)) = (w(X, tv),QF(X, tv + w(X, tv)} = 0. 

Therefore, for t 7^ 0, 

(v, F(X, tv + $(A, tv))) = t~ x {tv + w{X, tv),F(X, tv + w(A, tv))), 

and this is real by condition (|5|) of the theorem. Thus we can restrict 7 to a function 70 : R x R — > R. 

By a standard application of the Implicit Function Theorem to i _1 7o(A,i) = 0, in which (JTJ) - (|3j) are used (see for 
example 0), there is e > and a real-analytic function tfi\ : (— e, e) — * R such that <fi\(0) = Xq and if 70 (A, t) = 
with \t\ < e, then either t = or A = <fi\(t). Recalling that 7(A, e la t) — e la ^(X, t), we have shown that if 7(A, s) = 
and \s\ < e, then either s = or A = ^(|s|). 

We also note that by the symmetry, <fi\(—t) = <fi\(\t\) — 4>\(t), so <j)\ is an even real-analytic function, and 
therefore must in fact be a function solely of \t\ 2 . We therefore set 4>\{t) = (j)\(Vi), and so <p\ is real-analytic. 

We now define <f>^, : (— e, e) — >• R to be 

M*) = | t = 0j (40) 

is also real-analytic and satisfies s(j)^(\s\ 2 ) = w((j)\(\s\ 2 ), sv) for any s G C with \s\ 2 < e. 

Now we know that there is a neighbourhood of (Ao,0) in R x nullD^,F(Ao, 0) such that in that neighbourhood 
F{X,%b) = if and only if 7(A, s) — where Pip = sv. By taking a smaller neighbourhood if necessary, we have 
proven that F(X,ip) = in that neighbourhood if and only if either s = or A = 0a(|-s| 2 ). If s = 0, we have 
if) = sv + s<fi^(\s\ 2 ) = which gives the trivial solution. In the other case, ijj = sv + s(j>jf,(\s\ 2 ) and that completes the 
proof of the theorem. □ 

B The Operators L and M 

In this appendix we prove Theorems [6] and [5] The proofs below are standard. 

Proof of Theorem The fact that M is positive follows immediately from its definition. We note that its being 
strictly positive is the result of restricting its domain to elements having mean zero. □ 
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Proof of Theorem^ First, we note that L n is clearly a positive self-adjoint operator. To see that it has discrete 
spectrum, we first note that the inclusion H 2 L 2 is compact for bounded domains in R 2 with Lipschitz boundary 
(which certainly includes lattice cells). Then for any z in the resolvent set of L n , (L n — z) -1 : L 2 — > H 2 is bounded 
and therefore (L n — z)^ 1 : L 2 — » i 2 is compact. 

In fact we find the spectrum of L n explicitly. We introduce the harmonic oscillator creation and annihilation 
operators 

n tn 

L\ = d Xl T id X2 T + y^- (41) 
One can verify that these operators satisfy the following. 

1. = In. 

2. L n -n = -L n + L n _. 

3. L n - n = -L"lL\. 

As for the harmonic oscillator (see for example [U]), this gives the explicit information about <r(L) as stated in the 
theorem. 

For the dimension of the null space of L, we need the following lemma. 
Lemma 16. null(L" — n) = nullL" . 

Proof. If L™ip — 0, we immediately have (L n — n)ip = —LftLPip — 0. For the reverse inclusion we use the fact that 
||L"V>|| 2 = (LI^Llip) = (ip^L^ULtp) = (ip,(L n -n)ip). □ 

We can now prove the following. 

Proposition 17. nullL™ — n is given by 



nuU(L" - n) = { e ^ x ^ Xx+ix ^ ^ ^^VStaTfu+ira) | Cfc+?i = e ^^ e 2 «^ C)fc }. (42) 

k— — oo 

and therefore, in particular, dime null X™ = n. 

Proof. A simple calculation gives the following operator equation 
e fMV e -f|*| 2 =9xi + ^ 

This immediately proves that -0 € nullX" if and only if £ = e^l 2 ?/; satisfies c^i; + id X2 l; = 0. We now identify 
x G M 2 with z = x 1 + ix 2 G C and see that this means that £ is analytic. We therefore define the entire function O 
to be 

O(z) = e £ I — 

The quasiperiodicity of tp transfers to as follows. 

6(z + 7r) = e(z), 

Q{z + ttt) = e- 2mz e~ in7rTZ Q(z). 

To complete the proof, we now need to show that the space of the analytic functions which satsify these relations 
form a vector space of dimension n. It is easy to verify that the first relation ensures that have a absolutely 
convergent Fourier expansion of the form 



6(z)= J2 c ^ 



2kiz 

<^kC 
k— — co 

The second relation, on the other hand, leads to relation for the coefficients of the expansion. Namely, we have 

Ck+n = c m ^e 2k ^c k 

And that means such functions are determined solely by the values of Co, . . . , c„_i and therefore form an n-dimensional 
vector space. □ 

This completes the proof of Theorem [5] □ 
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C Fixing the Gauge 

We provide here an alternate proof of Proposition 02 largely based on ideas in [14]. We begin by defining the function 
B : R -> M to be 



B(C) = - f curl AfoQdt. 
r Jo 



It is clear that b = -^r J^ T2 B(Q d(. A calculation shows that B(( + rr 2 ) = B((). 
We now define P = (Pi, P 2 ) : R 2 -> K 2 to be 

a(ao = kc 2 - / 2 s(c)dc, 

Jo 

P 2 (x) = /" 1 curM(f,X2)de + ^^S(s2). 

A calculation shows that P is doubly-periodic with respect to C. 
We now define rf : M 2 -> R to be 



6 



r ] '(x) = -x 1 x 2 - AifatijdS- A 2 (x 1 ,0-P 2 (x 1 X)dC- 
1 Jo Jo 

rf satisfies 

Vr? = -A + A + P. 

Now let 77" be a doubly-periodic solution of the equation Arj" = — divP. Also let C = {C\, C 2 ) be given by 

where 17 is any fundamental cell, and set rj'" = G\X\ + C 2 x 2 . 

We claim that rj = rf + rf' + rf" is such that A + V?y satisfies (i) - (iii) of the proposition. We first note that 
A + V77 = A-A + A + P + Vt?" + C. By the above, A' = P + V77" + C is periodic. We also calculate that 
div A' = div P + A?/' = 0. Finally / A' = j P + Vry - C = 0. 

All that remains is to prove (iv). This will follow from a gauge transformation and translation of the state. We 
note that 

A Q (x + t) + A'(x + t) = A (x) + A'(x) + h -(^ ~* 2 

This means that A (x + t) + A'(x + t) = Aq(x) + A'(x) +Vg t (x) 1 where gt{x) = ^tAx + C t for some constant C t - To 
establish (iv), we need to have it so that Ct = for t = r, rr. First let I be such that r A I — —^f- and rr A I = —^f^- 
This I exists as it is the solution to the matrix equation 

r \ f h \ ( -f 



-i'T-1 rr :i J \ h ) \ % 



and the determinant of the matrix is just r 2 r 2 , which is non-zero because (r, 0) and rr form a basis of the lattice. Let 
C(x) = |i A x. A straight forward calculation then shows that e l ^^tp(x + I) satisfies (iv) and that A(x + 1) + VC(x) 
still satisfies (i) through (iii). This proves the proposition. 
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